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Abstract. Neutral hydrogen around high-z QSO and an optical depth t ~ 0.17 can be reconciled if reionization is more 
complex than a single transition at z ^ 6-8. Tracing its details could shed a new light on the first sources of radiation. Here we 
discuss how far such details can be inspected through planned experiments on CMB large-scale anisotropy and polarization, 
by simulating an actual data analysis. By considering a set of double reionization histories of Cen ( 2003 ^ type, a relevant 
class of models not yet considered by previous works, we confirm that large angle experiments rival high resolution ones in 
reconstructing the reionization history. We also confirm that reionization histories, studied with the prior of a single and sharp 
reionization, yield a biased t, showing that this bias is generic. We further find a monotonic trend in the bias for the models that 
we consider, and propose an explanation of the trend, as well as the overall bias. We also show that in long-lived experiments 
such a trend can be used to discriminate between single and double reionization patterns. 
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1. Introduction 

The first-year WMAP data release' detected a strong 
anisotropy-polarization cross-correlation at low / (Kogut et 
al. 2003). Previous analyses suggested that the Intergalactic 
Medium (IGM) had reionized at a redshift z - 6-8, and led one 
to expect an optical depth for Thomson-scattering r ^ 0.05 
(Miralda-Escude 2003 1. Observations of the Gunn-Peterson 
elfect in high-z QSO, requiring a fraction of neutral hydrogen 
at z ~ 6-7, agreed with this scheme (Djorgovski et al. 2001 
Becker et al. 2001 ; see however Malhotra & Rhoads 2004). On 
the contrary, the level of anisotropy-polarization correlation, in 
WMAP data, indicates r « 0.17 (Kogut et al. i2003l) and reion- 
ization at z ^ 16, assuming a single-step reionization model. If 
so, the reionization history could have been quite complex. In 
turn, data on CMB anisotropy and polarization could shed new 
light on the birth and evolution of primeval objects 

Different options on the nature of primeval objects have 
been considered. Ciardi et al. (2003") showed that metal-free 
stars in early galaxies may account for a depth up to r a; 0.15, 
although this seems an upper limit to such a picture (Ricotti 
& Ostriker "2004a). Various authors also considered an early 
pre-ionization due to black holes in small galaxies (Ricotti 8c 



Ostriker '2004b^ or miniquasars (Madau et al. '2004V or even 
the effects of sterile neutrino decay (Hansen & Haiman l2004> . 

Independently of the ionizing mechanism, however, com- 
paring WMAP and QSO data suggests that reionization is not 
achieved in a single, rapid step, but involves at least two differ- 
ent stages. Double reionization models were suggested, even 
before WMAR by Cen (2003 ), Wyithe et al. (2003 ), Sokasian 
et al. ( 2004J . Ricotti & Ostriker t2004 a) and others. In these 
models the ionization fraction attains a value » 1, at some 
high z, then partial recombination occurs, while a second reion- 
ization takes place at z = 6-8 and then x^. = 1 (for discussion 
purposes, we neglect Helium ionization). 

Haiman & Holder (2003 ) also considered another option, 
that the Universe partially reionized (up to Xp ~ 0.5-0.8) at a 
high Zr, to achieve a complete reionization at z - 6-8. A sim- 
ilar 2-step reionization was also schematically considered by 
Kaplinghat et al. 120031 for its effects on CMB. Still different 



reionization patterns were treated by BruscoU et al. ( I2002> . Hu 
& Holder (2003 ), Naselsky & Chiang ( l2004l and Colombo 
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I l2004t . confirming the interest in the relation between CMB 
data and early reionization, although often specifying no link 
with early object formation models. The present work is based 
on the dependence of the low-Z behavior of the Cf spectra on 
the whole reionization history, and aims to predict the actual 
detectability of r, as well as of the reionization redshift(s), the 
ionization rate(s), etc., with data on large angular scales. When 
dealing with full-sky small-angle experiments, the high num- 
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ber of pixels forces one to work in harmonic space which, under 
the assumption of a Gaussian signal, allows a fair information 
compression. For large-angle data, working directly in pixel 
space, as opposed to harmonic space, is numerically feasible 
and allows one to take into account a number of features that 
experiments can hardly avoid. 

In particular, the actual CMB signal can be recovered only 
on a portion of the celestial sphere. An example of such 
a limitation is the effect of Galactic contamination, making 
anisotropy data unreliable within ~ 20° of the galactic plane. 
On the other hand, first year WMAP data indicate that the main 
polarized foreground at frequencies up to 70 GHz is Galactic 
synchrotron (Bennett et al. 2003 ), whose polarization spectrum 
has a steep dependency on frequency. According to the syn- 
chrotron template by Bernardi et al. (2004), a lesser contam- 
ination is expected in polarization and the analysis can avoid 
Galactic cuts for such a signal. 

Combining temperature and polarization data covering dif- 
ferent sky areas is straightforward when working in pixel space, 
an important feature due to the different Galactic cuts possible 
for these two data sets. Taking them into account in the har- 
monic space requires either analytical approximations, which 
do not suite large-angle experiments well (see, e.g., Ng & Liu 
11999 ). or extensive Monte-Carlo simulations to calibrate suit- 
able window functions. 

In this work we perform a likelihood analysis of double 
reionization models of the Cen (2003| type. The reionization 
history is described by two parameters: besides r, which is re- 
garded as the reionization parameter more directly constrained 
by experiments, we consider z, , the redshift at which the IGM 
ionizes for the first time. A second reionization is then assumed 
to occur at z = 7. The duration of the first ionized period is 
fixed once by assuming that, between the two ionized eras, 
Xe = 1/3. This category of reionization histories is meant to ap- 
proach the pattern indicated by Cen ( 2003 ). We consider a grid 
of models spanning a large portion of the parameter space, and 
take into account temperature, polarization and temperature- 
polarization cross-correlation spectra. 

All of our analyses have been performed considering the 
features of the SPOrt^ experiment, including its sky coverage 
which avoids the Celestial polar caps. In addition, we have sim- 
ulated higher sensitivities in order to set experiment require- 
ments to significantly measure both relevant parameters, taking 
into account Cosmic Variance (CV). 

Cosmic reionization has its main impact on the harmonics 
linked to angular scales subtending the cosmological horizon at 
reionization. Accordingly, a knowledge of spectra above / ~ 30 
has a modest impact, so that the effect of using an angular res- 
olution ~ 5°-7° is also modest. The resolution of SPOrt (~ 7°) 
allows then a fair inspection of polarization features related to 
this epoch (patchiness effects, altering I ^ 1500 multipoles, are 
not considered here). A wide beam implies that beam smooth- 
ing must be taken into account in data analysis (on the contrary, 
the low-/ multipoles are almost free of beam effects in high- 
resolution experiments). Here we verify that an accurate treat- 
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ment implies no serious limitation in reconstructing the history 
of physical events, in a 7° experiment. 

Similar analyses have been performed in previous works 
using different reionization histories. Kaplinghat et al. (120031 
modeled a two-step reionization, while Holder et al. J2003> 
simulated reionization histories based on physical models of 
formation and evolution of the first ionizing sources. At vari- 
ance with the double reionization model by Cen ( I2003t and 
with models considered here, almost all models discussed in 
these works display a monotonic increase of the ionized frac- 
tion with time. 

All of these analyses coherently imply that, with levels of 
sensitivity similar to WMAP, information on reionization going 
beyond the integral t value is hardly obtainable. An important 
finding of Holder et al. (2003 ) is that t estimates resulting from 
fitting a sharp reionization history to models characterized by 
more complex reionization, are affected by a bias. Their nu- 
merical analysis shows that, at the WMAP sensitivity level, the 
estimate of t lays within 1-cr from the actual value. However, 
they show that the bias grows with sensitivity and the discrep- 
ancy between the actual and the estimated values can even go 
beyond ~ 10 cr for CV-limited experiments. 

This behavior is found also for the double-reionization pat- 
terns considered in this work. We also find that, for the noise 
range discussed here, the r estimate shifts to greater values as 
sensitivity improves. In particular, at sensitivity ~ WMAP, the 
optical depth tends to be underestimated, while a decrease of 
noise by a factor ~ 10 leads to an overestimated t, well outside 
3 cr, if a sharp reionization is assumed. We also provide an in- 
terpretation of why this trend occurs, and suggest that a test of 
the shape of the reionization pattern can be performed, in long 
term experiments, by comparing t estimates performed on data 
obtained in shorter or longer periods. 

The paper is organized as follows: in section|2lwe provide 
more details on the reionization histories that we considered 
and describe how the publically available code CMBFAST was 
modified to allow power spectrum computation for such histo- 
ries. In section |3] we describe the likelihood analysis. In sec- 
tion|4]we discuss the results of the analysis and in sectionlsjwe 
draw our conclusions. 

2. Reionization pattern and E-mode spectrum 

When the Universe reionizes, the CMB photon distribution has 
a relevant quadrupole term, greatly enhanced with respect to 
the last scattering epoch. This enables Thomson scattering po- 
larization to be preserved in the distribution of scattered pho- 
tons. Accordingly, the CMB polarization rate (CMBP) depends 
on the scattered fraction and, therefore, on t. The polariza- 
tion distribution on spherical harmonics, instead, depends on 
the evolution of the ionized fraction, so that the polarization 
spectrum bears an imprint of the reionization history (see, e.g., 
Zaldaii-iaganW KapHnghat et al. l2003l Holder et al. l2003l 
Naselsky & Chiang 2004). A quantitative evaluation of these 
effects can only be done numerically. Here we report results 
obtained by suitably modifying the publically available linear 
code CMBFAST (Seljak & Zaldarriaga .l996j . 
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Fig. 1. The parameter space allowed in our investigations in 
the plane r-vs-Zr- Dots show the position of fiducial models 
analyzed in this work. Fiducial models falling on the bottom 
edge have a single reionization. 

We select a flat ACDM cosmology with Hubble parameter 
(in units of 100 km/s/Mpc) h - 0.71 and density parameters 
n„h^ = 0.148, Qfc/i^ = 0.024. We then take a grid of points 
in the t-z,- parameter plane, spanning the intervals 0.07 < r < 
0.30 and 10 < z, < 39 (see FigureU). 

Double reionization allows the same value of t with dif- 
ferent Zr- The minimal allowed Zr for such a r value is obtained 
for a single sharp reionization occurring at a suitable z; this sets 
the lower side of the shaded area in Figure [2 Then, the low- 
ionization period, lasting from a redshift zi to zi =7, reduces 
to zero. For any greater Zr, 

T^T2+ 0-TCt2 [- ) --(- ) -- (1) 

Here o-j- is the Thomson cross-section, nip is the average 
baryon mass, phizi) the baryon matter density at redshift zi', 
t2 is the time corresponding to Z2, T2 is the optical depth due 
to full reionization since f2. In eq. (Q} a matter dominated ex- 
pansion is assumed, since Zr to Z2 - 7, yielding an error of a 
few percent. Given t and the (1 +Zr)/(1 +22) ratio, eq. Q fixes 
the ratio (1 + zi)/(l + Z2)- There is however a top value for Zr, 
achieved when Zr - Z\; this sets the upper side of the shaded 
area in Figuren(obtained through exact numerical integration). 
Then, the first Xe - 1 period, lasting from Zr to z\ , reduces to 
zero. 

A first reionization occurring at a redshift Zr » 30 is un- 
likely under most models of ionizing sources (see, e.g., Haiman 
& Holder ISOOSl l: therefore, not all points falling within the 
shaded area of Fig.^bear the same physical relevance, in spite 
of being compatible with the parametrization of reionization 
discussed here. 

When modifying the CMBFAST linear code, we are very 
careful in treating on ionization transients. Analytical and nu- 
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Fig. 2. E-mode angular power spectra (top panel) and ioniza- 
tion history (bottom panel) for three different models consid- 
ered in our analysis with t = 0.17. Solid line represent a ion- 
ization profile with Zr - 18, long-dashed line has Zr - 22 and 
short-dashed line Zr - 26. Different ionization levels are con- 
nected by a fast but smooth transition in order to guarantee sta- 
bility in numerical integration. 



merical approximations at each shift are the same as used in 
CMBFAST when dealing with a single reionization event. In 
particular, to avoid instabilities in numerical integration, due to 
fast variations in x^., we adopt grid steps Sj - .01 and - 1. 

In Figure|2] we show the varying ionization rates and the re- 
sulting CMBP /i-mode spectra for a set of models with equal t 
but different z, . Conversely, in Figure|3l we show models where 
the first reionization redshift z, is kept constant and t is vari- 
able. The models of Figure|2lie on a vertical line of Figure^] 
those of Figure |3] lie on a horizontal line in the same figure. 
Power spectra are clearly sensitive to the ionization history and 
potentially encode information on it. Notice that varying t has 
eff'ects extending to greater I, while varying Zr produces effects 
limited to I $ 50. In the next sections we discuss how far these 
variations are detectable at different levels of instrumental sen- 
sitivity. 

3. Likelihood Analysis 

Large angular scales are significantly affected by CV and a 
given model can yield significantly different skies. To evaluate 
how far large angle CMB experiments can recover the reion- 
ization history, we adopt a Monte Carlo approach. 

The basic outline of our approach is as follows. We select 
a fiducial cosmological model M. and generate 5000 sky maps, 
with resolution and beam smoothing similar to those of actual 
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Fig. 3. Cf power spectra and ionization histories for three mod- 
els with Zr - 25 and t = 0.17 (solid line), t - 0.21 (long- 
dashed line), T - 0.25 (short-dashed line). 



experiments. To each CMB map we add a noise map; the re- 
sulting maps for the anisotropy T and the Q, U Stokes param- 
eters are then used for likelihood evaluation. The probability 
distribution function (p.d.f.) for the parameters characterizing 
the model Ai is then simply the distribution of the best-fit pa- 
rameters in the likelihood analysis of such maps. 

In this work we deal with polarization data similar to those 
expected from the SPOrt experiment, however allowing for a 
higher sensitivity, while temperature data are supposed to come 
from an experiment similar to WMAP. Simulated maps are 
generated using the HEALPix^ package. Both T and Q, U data 
are smoothed with a 7° FWHM Gaussian filter and we chose 
a HEALPix resolution Nside - 16 (corresponding to a pixel 
width of ~ 3.5°). This approach results in a rather low number 
of pixels, while parameter extraction is sensitive only to Z ^ 40 
multipoles. Each simulated map represents a random realiza- 
tion of the process Ai + N. The noise model N is assumed to 
be uniform and white, and is then fully defined by the rms noise 
values cTj- and crp for T and Q, U pixels. 

As explained in the Introduction, we remove the region 
with Galactic latitude \b\ < 20° from T maps, where Galactic 
contamination strongly surpasses the CMB signal. On the other 
hand, the synchrotron template developed by Bernardi et al. 
(2004) shows that at 90 GHz the synchrotron polarized emis- 
sion in the Galactic Plane is at most comparable with the CMB 
polarization signal for a cosmological model with optical depth 
T = 0.17. Taking into account that foreground removing tech- 
niques allow one to lower the foreground contamination by a 

^ http://www.eso.org/science/healpix/ 



significant amount (by a factor ~ 5-10), the residual contam- 
ination is likely to affect the CMB polarization analysis to a 
negligible extent, and we choose not to remove the Galactic 
Plane in the present analysis of (Q, U) maps. Still, part of the 
Galactic Plane is cut out by excluding declinations \6\ > 51.6°, 
which SPOrt is unable to inspect. (In addition, we notice that 
polarized emission by dust grains is expected to lie safely be- 
low synchrotron, at least up to Z ^ 50; see Fabbri 120041 for a 
more detailed discussion). 

Simulated data are ordered into a vector x = {T(i - 
\,...,NT),Q{i= \,....,Np),U{i= \, ....,Np)], Nt Npheing 
respectively the number of anisotropy and polarization pixels. 
The likelihood of a model M is then given by a multivariate 
Gaussian: 

X(A<W=^5;^^exp(-i,'C-4); (2, 

where the correlation matrix reads: 

C,.y = <xfx,.) = S,-,-+N,.; (3) 

(brackets mean ensemble average). Its elements represent the 
expected correlation between elements of the vector x. For ex- 
ample, the expected correlation between T signals in the i and 
j pixels, at an angular distance reads 

iTiTj) = '^'V)'B? + i'^rfSif, (4) 

here is the anisotropy spectrum of the model, f ;(cos 6) are 
Legendre polynomials, while the coefficients Z?^ account for 
pixelization and beam smoothing. For a FWHM of 7°, the cor- 
rection is relevant even for the lowest harmonics. Similar 
expressions, taking into account the tensor nature of polariza- 
tion, hold for correlations involving Q and U data (see, e.g., 
Zaldarriaga 1998 Ng & Liu 1999 1. The effects of sky cuts, 
here, are directly set by the /, j index domains. This is a clear 
advantage of working in the coordinate space. (As already out- 
lined, in the harmonic space, fsky corrections are unsafe for 
large angle experiments or when T and Q, U data cover dif- 
ferent portions of the sky.) 

For each sky realization, we seek the values of the reion- 
ization parameters which maximize the likelihood function X. 
Repeating the operation for a set of realizations allows us to 
study the resulting distributions. 

The fiducial models considered in this work are shown in 
Figure ^ Each model was tested at 3 levels of sensitivity for 
the polarization measures: crp = 1.50, 0.45, 0.15yL(K for ~ 7° 
pixels, corresponding to ~ 10, 3, 1 //K-degree. The first value 
is the expected sensitivity for the 90GHz SPOrt channel. The 
pixel noise for T data is set at the reference value crp = IfiK. 
For such cTp, uncertainties on measurement of T-multipoles 
are dominated by the CV up to Z ^ 500 while, with our 
beamwidth and pixelization, only the first ^ 40-50 multipoles 
matter Thus, a ctj reduction would not yield a significant gain. 
On the contrary, as polarization is ~ 100 times smaller than 
anisotropy, sensitivity increases like those considered here im- 
prove the ratio between CV and noise variance, mainly in the 
harmonic range we are exploring. How the different crp values 
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Table 1. Fiducial models: r is the total optical depth, Zr is the 
redshift of the first reionization. Models with the superscript s 
have a sharp reionization history, other models assume that the 
Universe is completely reionized at z < 7, while between these 
two reionization periods drops to 1/3. 
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Fig. 4. Distribution of r for models £), at different noise sensi- 
tivity (left to right, crp = 1 .50, 0.45, 0.15 yuK). SoHd lines show 
results obtained marginalizing over Zr, dotted lines refer to re- 
sults of checking the same realizations against single reioniza- 
tion models. At high sensitivity, the latter prior give rise to a 
noticeable bias. 



interfere with polarization multipoles must be considered in de- 
tail, to understand our results, and will be discussed below. 

The distribution of the best fit parameters, among the 5000 
CMB realizations of each fiducial model, at the three noise lev- 
els, provides a frequentist estimate of the probability density 
function (p.d.f.) of model parameters. To ensure that CV did not 
introduce spurious effects in the comparisons of different mod- 
els, we used the same set of random seeds and of noise maps 
for all fiducial models. Therefore: (i) Differences between sets 
of sky maps with the same noise level are due only to changes 
in the CMB power spectra between the fiducial models, (ii) 
Differences between sets of maps at different sensitivities (but 
corresponding to the same model) are only due to variations in 
the noise. 

4. Results 

All fiducial models were analyzed first assuming a single reion- 
ization, then considering reionization histories of the kind de- 
scribed in Sec. |2] We compare the probability distributions on 
T, to test the bias induced by the prior of sharp reionization. 

In Fig. 13 we compare the t distributions obtained under 
the prior of single reionization (dotted lines) with those ob- 
tained after marginalizing over Zr the joint distribution on r 
and Zr (solid lines). The three plots account for three polar- 
ization noise levels. At each level, no significant difference be- 
tween the width of the two distributions is apparent. However, 
while all marginalized distributions peak at the true optical 
depth T - 0.17, the single reionization distributions show a 
bias, depending on crp, similarly to Holder et al. 12003) . Let 



Fig. 5. E-mode angular power spectra for models C (solid 
lines), D (long-dashed lines), fi (short-dashed) and two sharp 
reionization models with t = 0. 15, 0. 19 (dot-dashed and dotted 
lines respectively). At given value of the optical depth, spec- 
tra of models with double reionization fall below the APS of 
a sharp reionization with same r for / < 8-10 and above for 
10 < Z < 40. The long-short dashed lines, instead, show the 
noise spectra for different pixel noises on polarization (top to 
bottom (Tp - 1.50, 0.45, 0.15;uK), corrected for pixelization 
and beam smoothing. 



us however outline a further trend: if crp = 1 .50 f^K, the single 
reionization prior leads to underestimating t, although within 
the (still wide) statistical error; as crp decreases, the peak r- 
value increases and, at the top sensitivity considered, t is over- 
estimated by more than two standard deviations. 

This happens because the range of Cf significantly affect- 
ing the estimates is strongly related to crp, for the sensitivity 
levels we are considering. In fact, t accounts for the number of 
scattered CMB photons, which can affect the total amount of 
large scale polarization (Zaldarriaga 1997p. Therefore, the E- 
spectra, for models with different r but equal Zr, differ mainly in 
the height of the first reionization peak. But a variation of z, - 
at fixed t - also shifts the peak in /. Fig.|5]shows the E-spectia 
for the models C, O, &, plus two single reionization models, 
with T - 0.15 and t - 0.19. For models with t - 0.17, when 
Zr increases, the first reionization peak moves to greater I and 
its height slightly decreases: then / % 8 harmonics of models 
£) and & approach those for single reionization with lower t; 
the same models, in the range 8 ^ / $ 20, resemble a single 
reionization with greater t. 

When crp - \ .50 //K (top noise level considered), reioniza- 
tion parameters are mostly determined by the first 6-7 multi- 
poles, the only ones clearly above noise (see Fig.|5}. Trying to 
fit a single reionization to a double reionization history there- 
fore leads to underestimating t. As the noise decreases and 
multipoles in the range 10 ;S Z :$ 20 acquire a greater weight, 
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a double reionization history can be misinterpreted as single 
reionization with a greater t. 

The shift of the reionization peak to greater I is typical of 
models with complex reionization histories. Hence, this low-cr 
high-T effect, when a single reionization is assumed, should be 
independent of reionization details, for the range of sensitivi- 
ties considered here. For higher sensitivities, when multipoles 
20 :$ / 40 affect the r determination, the /s-spectrum exhibits 
further peaks, whose details are related to model features, and 
the trend will depend on the set of models considered. 

In general, the bias is stronger for greater r and, at fixed 
T, it increases with z,-. At noise levels accessible to current ex- 
periments' statistical errors exceed the bias, and r estimates 
are safe, although systematically smaller than the real values. 
In higher sensitivity experiments, wrong priors can cause a t 
overestimate, exceeding several standard deviations. An oppo- 
site behavior is found when single reionization models are an- 
alyzed assuming a two-reionization history. In this case, how- 
ever, the bias is not so severe. 

The low-cr high-T effect can be used to test the reionization 
pattern, by comparing early and late outputs, in a long lived ex- 
periment. Detecting a clear trend requires a sensitivity increase 
by at least a factor S 3. As we are considering low multipoles, 
early outputs must cover the whole sky with a reasonable sen- 
sitivity. Assuming that a full sky coverage requires ~ 6 months, 
the total experimental lifetime should be ~ 5-6 years. The re- 
quired increase is unlikely within the planned lifetime of cur- 
rent experiments, but could be an option for future polarization 
missions. 

We performed some more detailed tests of this point, (i) For 
each realization we labeled t\^^, t'^^^^t^^^ (rf^Q, t^45T^i5) 
the optical depth estimated at cTp = 1.5, 0.45, 0.15 yuK, under a 
single (double) reionization prior (ii) We selected the realiza- 
tions where the r estimates, under a single reionization prior, 
increased at both noise reductions (tJ^q < "^"045 < '''0 15)- 
Among them we kept those for which the shift Tj, - r j 
exceeded twice t^^^ - t^^q- This residual fraction averages 
~ 65%; although varying from ~ 20% (model S) to more than 
90% (model 'H). A progressive shift to higher values in single- 
reionization t estimates, when sensitivity increases, therefore 
can be considered as a hint that a more refined description of 
reionization is needed. 

An additional test, for long lived experiments, is provided 
by dividing the global dataset into subsets corresponding to 
shorter observation periods, and comparing the r value esti- 
mated by the analysis of the full dataset with the weighted av- 
erage of the estimates in the smaller subsets. For reference pur- 
poses, we considered a SPOrt-like experiment lasting 5 years, 
and achieving a polarization pixel sensitivity crp = 0.40yuK 
(Nside - 16). Neglecting correlated noise effects and assum- 
ing that a full sky coverage requires 6 months of observations, 
the global data can be divided into 10 smaller subsets, charac- 
terized by cTp ~ 1 .25yuK. We then compared the global r esti- 
mate with the average of the estimates in the smaller subsets. 
Repeating this analysis for 100 random realizations of model 
D, we found global and averaged estimates to be consistent 
with each other and with the actual r, when correct priors are 
made (top panel in Fig.|6ll. On the contrary, when sharp reion- 
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Fig. 6. Top. Open circles: t from a WMAP+ SPOrt-like exper- 
iment, with cTp - 0.40/iK in 5 years, for 100 realizations of 
model £), and with a fair double reionization prior. Crosses: 
average t obtained by dividing the whole data set into 10 sub- 
sets (6 months observation time). Error bars yield the 1-cr c.l. 
Bottom. The same as Top panel but assuming a sharp reioniza- 
tion prior. Global (averaged) estimates are shifted up(down)- 
ward in respect to the actual t value. 



ization is assumed, averaged estimates are systematically lower 
than the correct value, while global measurements overestimate 
it (bottom panel in Fig.|6j. Then, in Fig.0 we plot the distribu- 
tion of the ratio between the full 5-year estimate and the aver- 
aged 6-month estimates. Assuming a correct prior, such ratio 
averages T^yrliTbm) - 0.99 ± 0.07, while for sharp reionization 

T5yrl{T(„n) = 1.20 ±0.10. 

We now adopt a double reionization prior, and discuss the 
precision by which t and Zr can be recovered. Marginalizing the 
joint probability distribution in the T-Zr plane on either param- 
eter provides the ID probability density for the other parame- 
ter. The variances of these distributions, averaged over fiducial 
models, tell us how far the models can be discriminated at each 
sensitivity; results are displayed in Table|5] two equal-r models 
can be distinguished if their Zr are farther apart than ~twice the 
value shown. These uncertainties are close to those obtained in 
the analysis of the performance of high resolution experiments 
with equivalent sensitivities (see, e.g., Kaplinghat et al. l2003> . 
Thus, low angular resolution is not a serious impediment in 
reionization studies. 

Within the allowed region of parameter space (see Fig. ^ 
and discussion), Table|2]implies that r is better fixed than Zr, if 
the priors on the reionization history are correct. This confirms 
that also for double reionization, optical depth is the most rel- 
evant parameter. At a fixed noise level, models with higher r 
generally allow for a better estimation of both parameters, due 
to their higher signal-to-noise (S/N) ratio; among models with 
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Fig. 9. Same as Fig.|8] but for t after marginalization over z,-- 




Fig. 7. Distribution of the ratio between global and averaged r 
estimates in 100 realizations of model D, assuming either dou- 
ble realizations (solid hne) or sharp reionization (dotted hne) pig. 10. Same as Fig. |9l but for models £ (solid lines) and g 



priors. 



(dotted lines). 



Table 2. Average statistical errors for the recovered parameters 
as a function of polarization sensitivity crp (fiK). 
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Fig. 8. Distribution of maximum likelihood Zr, after marginal- 
izing over T, for 5000 realizations of models C (dashed lines), 
D (solid lines) and & (dotted lines). From left to right, panels 
refer to crp - 1.50, 0.45, 0.15^K, respectively. 



the same optical depth those with lower Zr have slightly sharper 
redshift distributions, because for very early first reionization, 
the differences between spectra with distinct Zr are less pro- 
nounced. 

In more detail, Fig.|S]shows the distribution of z,- for models 
with T - 0.17 (see caption for details). For crp - 1.50/iK (left 
panel), the probability distributions overlap significantly, and z,- 
cannot be recovered, while for crp - 0.15/iK (right panel) they 
are almost completely distinguished. The middle panel displays 




Fig. 11. Same as Fig.QOI but for Zr- 



an intermediate situation where models with single reionization 
can be distinguished from double reionization models with z,. in 
the upper half of the allowed range (i.e. Zr <: 25 for r - 0.17). 
Moreover, distributions on Zr, in double reionization models, 
are wider by a factor ~ 2 than in single reionization models. 
Figure 121 shows that the probability distributions on r, in mod- 
els D and £, are indistinguishable at all sensitivities; thus, dif- 
ferences in Zr do not greatly affect the ability to recover t, if 
correct priors are used. Model C, instead, displays a bias, as 
discussed above. In this case, t is progressively underestimated 
for decreasing crp, as a single reionization model is analyzed 
assuming a double reionization prior. 

In Fig. [To] we plot the distributions on r for models with 
equal Zr but different r (£ and Q). Even at the highest noise, a 
difference can be seen; for crp - 0.45 /iK the two distributions 
are clearly separated. In both cases, the correct t is recovered. 
Figure instead, shows results for Zr- The probability distri- 
butions of both models exhibit a significant overlap and cor- 
rectly peak at the true value of the reionization redshift. Model 
Q, however, is characterized by a sharper and better defined 
distribution, due to its higher S/N ratio. 
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We can summarize our findings as follows, (i) Constraints 
on the reionization history from large-angle experiments are 
similar to those by high-resolution experiments with equiva- 
lent sensitivity. For instance, a polarization noise crp - 1 .50 ijK 
allows one to constrain the total optical depth with an accuracy 
of ~ 20% for T = 0.17. Increasing sensitivity by an order of 
magnitude allows for a simultaneous detection of t and Zr, with 
an accuracy ^ 5%. (ii) A wrong prior on the reionization his- 
tory causes a bias in the r estimate. For crp - 1 .50 fiK the statis- 
tical uncertainty exceeds the bias, while for crp = 0.15 fiK the 
bias may exceed 3 standard deviations. This point, too, agrees 
with previous works studying the capabilities of high resolution 
CMB polarization measurements. Thus high resolution, by it- 
self, does not provide more reliable estimates, (iii) At current 
noise levels, r is underestimated, while, for increasing sensitiv- 
ities, it tends to be overestimated. This is likely due to the first 
reionization peak moving to higher /'s in models with complex 
reionization, with respect to sharp reionization models with the 
same t. 

(iv) Finally, we suggest an observational approach to dis- 
criminate between single and double reionization histories, 
based on such low-cr high-r effect. For illustration, we con- 
sidered a mission lasting 5 years, allowing a pixel sensitivity 
crp ~ 0.4 yuK. We assumed that full sky coverage requires 6 
months so that the full dataset can be divided into 10 smaller 
subsets, with a sensitivity worse by a factor ~ 3. Assuming 
no systematic effects, the 10 subsets provide as many indepen- 
dent T estimates. When the correct priors on reionization are 
assumed, the global estimate is consistent with the average of 
the estimates from the smaller subsets; on the contrary, assum- 
ing a sharp reionization prior, the 5-year t value systematically 
exceeds both the correct value and the average of the 6-month 
estimates. Accordingly, in future low-noise experiments, it can 
be significant to compare whole-run t estimates, with full sen- 
sitivity, with averages among short run estimates. Similar error 
bars are expected, but a shift of the estimated t would indicate 
that the assumed reionization pattern could be a source of bias. 



5. Conclusions 

In this work we discussed a class of physically motivated dou- 
ble reionization models, characterized by two parameters: the 
total optical depth t, and the first reionization redshift Zr- We 
determined at which sensitivity level their features can be re- 
covered by large angle CMBP measurements. We find that 
wrong priors on the history of reionization can lead to a bias 
in T estimates, in agreement with findings for high resolution 
experiments (Holder et al. l2003> . At the WMAP or SPOrt noise 
level, the bias is well within statistical errors; at higher sensi- 
tivities the T estimate can lie several standard deviations from 
its true value. Holder et al. (2003 1 argue that fitting a two-step 
reionization to the models they considered allows one to signif- 
icantly reduce the bias. In addition we find that, for the class of 
models considered here, the biased r estimates exhibit a char- 
acteristic dependence on experimental sensitivity. Testing this 
low-cr high-T effect in actual experiments provides an indica- 
tion of the fairness of priors. 



Within the context of double reionization models, a pixel 
noise crp = 1.50/iK for 7° pixels, allows one to constrain t 
with an accuracy ~ 20%; an increase in sensitivity by a factor 
~ 3 enables us to distinguish single reionization models from 
models with an early reionization (i.e. z,- ^ 25 for t - 0.17) 
followed by a partial recombination period. Models of ionizing 
sources rarely yield Zr ^ 30. Firm measurements of both t and 
Zr, with precisions ~ 5%, require a sensitivity increase by an 
order of magnitude. Greater spatial resolution, instead, is not 
very relevant, as most information on reionization is carried by 
the first 30^0 multipoles of the is-spectrum. 

Measures of Cf at large angular scales are an important 
probe of the evolution of ionizing sources at redshifts 10 ^ z :^ 
30. However, as finer data become available, the situation be- 
comes more and more risky. Unless an accurate class of models 
is fitted to data, parameter misestimates can occur: apparent er- 
rors can seem small, while true values lie well off the 3 cr error 
interval. A test of the reliability of the overall t estimate can 
be however performed by comparing estimates from the whole 
data set with those from smaller representative subsets. A clear 
detection of this low-cr high-r effect requires an increase in 
sensitivity by a factor ~ 3. Assuming that full sky coverage re- 
quires 6 months, the effect can be evident in experiments last- 
ing ~ 5-6 years. 
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